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We study the filling-controlled metal-insulator transition in the two-dimensional Hub- 
bard model near half-filling with the use of zero temperature quantum Monte Carlo 
methods. In the metallic phase, the compressibility behaves as k oc |/i — ^ ic |-°- 58±0 - 08 
where \x c is the critical chemical potential. In the insulating phase, the localization 
length follows oc \\x — ^ c \~ Ul with v\ = 0.26 ± 0.05. Under the assumption of hyper- 
scaling, the compressibility data leads to a correlation length exponent v K = 0.21 ±0.04. 
Our results show that the exponents v R and v\ agree within statistical uncertainty. This 
confirms the assumption of hyperscaling with correlation length exponent v = 1/4 and 
dynamical exponent z = 4. In contrast the metal-insulator transition in the generic 
band insulators in all dimensions as well as in the one-dimensional Hubbard model 
satisfy the hyperscaling assumption with exponents u = 1/2 and z = 2. 
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Physics of the metal-insulator (M-I) transition has been studied over decades from the 
point of view of quantum critical phenomena and strong correlation effects. Recently, an 
approach to the M-I transition using a hyperscaling hypothesis has been developed by one of 
the authors (M.I).B' Let us consider the filling-controlled M-I transition in the ground state, 
which is controlled by the shift of the chemical potential 

A = fi-fi c . (1) 

Here fj, is the chemical potential while /i c is its critical value at the M-I transition. If we 
assume the single parameter scaling hypothesis, the correlation length of the system £ is 
described as 

e~iAr, (2) 

and the singular part of the free energy is given byi'& 

/ S (A) ~ \A\^ d+z \ (3) 

Here, z is the dynamical exponent, while d is the dimensionality of the system. In the 
metallic phase, the doping concentration 5 = 1 — n near the phase transition is given by 

so that from eq. @ we have 

Here we used the relation vz — 1 which is obtained from the generalized Josephson relation. 

In the case of the Hubbard model on a square lattice near half-filling, it has been shown 
by two of the authors (N.F. and M.I.)Bfl (hereafter NF-MI) that 5 satisfies 

\A\oc5 2 . (6) 

This result was obtained by computing the chemical potential, 

dE G 

at U/t = 4 on 4 x 4 to 12 x 12 lattices, in the standard notations. Here Eq denotes the 
ground state energy. The calculations were carried out with the zero-temperature quantum 
Monte Carlo algorithm.SB* This result implies that the compressibility diverges as 

K =^~ ' A '- 1/2 - < 8 > 

when we approach the M-I transition from the metallic side. However, in NF-MI, an esti- 
mate of the statistical error of the exponent was not done, mainly due to large statistical 
uncertainty in location of the critical point, /i c . 
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Recently, two of the authors (F.F.A. and M.I.)§B (hereafter FFA-MI) have computed with 
QMC methods, zero temperature imaginary time displaced Green function in the insulating 
phase. From this data they obtain an accurate estimate of the critical chemical potential at 
U/t = 4: jj, c — 0.67 ± 0.05 in units of the hopping matrix element. This result is obtained 
by extrapolating to the thermodynamic limit from data on 4 x 4 to 16 x 16 lattice. 

Combining this value with the chemical potential data of NF-MI in the metallic phase 
leads to: 

|A| oc 5 2 - 38±0A3 (9) 

(see Fig. 1). From this result and the assumption of hyperscaling we obtain a value of the 
correlation length exponent 

v K = 0.21 ±0.04. (10) 




Fig. 1. Log-log plot of \[i — fi c \/t verses S in the metallic phase. Here, fi c /t — 0.67 ± 0.05. The line is the 
least-squares fit to the data. 
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In order to directly determine the correlation length exponent, FFA- have computed 
the Green function G(r,u = fx) for values of the chemical potential within the charge gap. 
For those values of the chemical potential, G(r, to = //) oc e~' r '/^. The localization length, £/, 
diverges as the M-I transition is approached from the insulating side. We obtain numerically: 

6 oc |A|-°' 26±a05 (11) 

(see Fig. 2). From this result, we obtain an estimate of the correlation length exponent: 

^ = 0.26 ± 0.05. (12) 
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Fig. 2. Log-log plot of \n — fJ, c \/t verses in the insulating phase. Here, fi c /t = 0.67 ± 0.05. The line is 
the least-squares fit to the data. 



If hyperscaling is valid, one expects the relation: 

v K = v h (13) 

That the exponents should be equal on either side of the critical point is obtained through 
the scaling properties, and from the fact that the insulating and metallic phases may be 
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connected smoothly around the critical point A = 0.I& The main result of this paper is 
that eq. (|13|) is satisfied within statistical uncertainty, thus putting on a firm numerical basis 
the assumption of hyperscaling. The M-I transition in 2D Hubbard model belongs to the 
universality class characterized by v — 1/4 and z — 4. This stands in contrast to the generic 
band M-I transition in all dimensions as well as to the Mott transition in the ID Hubbard 
model, which both belong to the universality class v = 1/2 and z = 2. The M-I transition in 
2D Hubbard model belongs to the universality class characterized by v = 1/4 and z = 4.!) 
This stands in contrast to the generic band M-I transition in all dimensions as well as to 
the Mott transition in the ID Hubbard model, which both belong to the universality class 
v = 1/2 and z = 2.0) 
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